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MOTIVATION 

Approximation schemes have played an important role 
in the development of General Relativity (GR) since its 
formulation back in 1915, cf. [1, 2, 3, 4, 5, 6, 7, 8, 9] 
for some early approaches. Several schemes have been 
devised to quantify relativistic gravitational effects in the 
context of isolated gravitating systems, thereby permit- 
ting a confrontation of the at that time new gravitational 
theory with experiments. Even many modern applica- 
tions, which take into account general relativistic effects, 
rely on approximative methods [10, 11]. 

Nowadays it is common practice to assess the success 
of new gravity theories by formulating them in a stan- 
dardized parametrized language [12], which allows for 
a rapid comparison with the outcome of several exper- 
iments. Despite the ample use of such a scheme in the 
context of isolated gravitating systems no such formula- 
tion exists in a cosmological context. The lack of such 
a scheme is rather surprising since, on a first guess, one 
would expect that cosmological scales demand for a rel- 
ativistic treatment and, therefore, for an approximation 
which allows us to classify different cosmological mod- 
els in a systematic way. 

In the following we outline the program which has to 
be carried out in order to formulate such a framework in 
cosmology and report on recent progresses in formulat- 
ing the first order post-Newtonian (1PN) cosmological 
hydrodynamic equations for Einstein's theory. 

APPROXIMATION SCHEMES 

Although there exists a whole host of different approx- 
imation schemes in relativistic gravity we are going to 
concentrate only on two in the following. At this point 
we also have to remind the reader that there is no com- 
monly accepted definition of these schemes, and, there- 
fore, the usage of the names for them is far from unequiv- 



ocal. In order to fix our usage of the terms we briefly in- 
troduce the so-called post-Minkowskian (PM) and post- 
Newtonian (PN) approximation. 

In the post-Minkowskian or weak-field / fast-motion 
approximation one extends the gravitational potential, 
i.e. the metric in Einstein's theory, around the flat 
Minkowski metric 1 g^ v — rj^ v + ey^ + .... In many ap- 
plications one then tries to rewrite the field equations of 
GR in a form which closely resembles the form of the in- 
homogeneous wave equation, cf. [13, 14, 15, 16, 17, 18]. 
The solution of the field equations then being constructed 
by the use of a retarded Green's function. No other matter 
variables are expanded in the PM approximation. 

The strategy in the post-Newtonian or weak-field / 
slow motion approximation is different. Here one devel- 
ops the metric, starting from the Newtonian limit of Ein- 

o 

stein's theory denoted by g^ v , into a series of inverse 

o t 

powers of the speed of light, i.e. g^ v =g^ v +c _1 g^ v 
+ .... In addition the four velocity, which enters the defi- 
nition of the energy-momentum tensor, is developed into 

o 1 

a series u^ 1 =u^ +c _1 + . . . in inverse powers of c. The 
strategy in solving the field equations in the PN approx- 
imation also differs from the one employed in the PM 
approximation. The goal in the PN approximation is to 
recover a set of equations which resemble the form of the 
field equation in Newtonian gravity, i.e. Poisson's equa- 
tion. Since GR has more structure than Newton's theory, 
one usually has to deal with a coupled set of Poisson like 
equations in the PN approximation. 

Different ways of labeling the post-Newtonian order, 
i.e. counting the powers of c up to which the metric 
components need to be determined, can be found in 
the literature. Here we stick to the counting scheme 
originally introduced by Chandrasekhar and coworkers 
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in a pioneering series of works on the PN approximation 
for isolated gravitating systems, cf. [19, 20, 21, 22, 23]. 



COSMOLOGY 

Let us now come to the situation in cosmology. As de- 
picted in figure 1 one currently applies different ap- 
proximations on varying cosmological scales. From a 
more general standpoint the approximations used in cos- 
mology can be traced back to the schemes described 
in the last section. The theory of cosmological pertur- 
bations, e.g., could also be viewed as a modified ver- 
sion of the post-Minkowskian approximation. The ma- 
jor difference is that one now starts from the cosmologi- 
cal Friedman-Lemaitre-Robertson-Walker spacetime in- 
stead of the Minkowski spacetime. 

It is interesting that there has been made ample use of 
PM like approximation techniques in cosmology. In con- 
trast an approximative scheme in the spirit of the previ- 
ously described post-Newtonian approximation has not 
attracted much attention during the last years [24, 25]. 
This is even more surprising if one thinks of impor- 
tant cosmological tests like the formation of structure, 
in which the physical assumptions of weak gravitational 
fields and slow motions find a direct application. In the 
next section we briefly comment on the form of the equa- 
tions of motion in the cosmological PN approximation. 

Finally, let us note that from figure 1 it also becomes 
immediately clear that in cosmology we need a unifica- 
tion of both schemes, i.e. PM and PN, into one approxi- 
mative framework which covers all scales. 



Post-Newtonian cosmology 

Recently we started the post-Newtonian part of the 
program in [26, 27]. In order to provide a basis for future 
numerical structure formation simulations we derived the 
cosmological equations of hydrodynamics up to the first 
post-Newtonian order (1PN). They read 2 : 
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The dot "•" denotes a partial derivative with respect to (. 



Here the quantities U,V, <5,V>, and ho a denote gravi- 
tational potentials from the different metric components. 
For details see [27]. As usual, the cosmological scale fac- 
tor is denoted by a and the quantities p , II , p , v a represent 
the matter variables of the system, i.e. energy-density, in- 
ternal energy-density, pressure, and three velocity of the 
fluid elements. 

The equations (1) and (2) replace the continuity equa- 
tion and Euler equations that are usually used in New- 
tonian studies. The 1PN equations should be used in fu- 
ture simulations to describe the first order corrections of 
General Relativity to cosmological hydrodynamics, for 
example in codes which study the formation of structure 
in the non-linear clustering regime. Although we expect 
that the post-Newtonian c~ 2 corrections to the equations 
of motion are small in many situations, we cannot pre- 
clude that there are secular effects in the structure forma- 
tion process. 

In other words, one still has to show, by performing 
a numerical simulation on the basis of the higher than 
Newtonian order equations, that the currently used New- 
tonian description is appropriate. 

The set of field equations to be solved in connection 
with the equations of motion (1) and (2) are given in [27]. 
The equations therein are kept in a form which closely 
resembles Chandrasekhar's approach [19] in the non- 
cosmological case. In [26] we also consider a mixture 
of post-Newtonian and perturbative methods and provide 
the final set of equations for several gauge conditions. 
At this point we would like to stress that, depending 
on the needs of the people who perform simulations, 
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FIGURE 1. The current cosmological tests make use of different approximation schemes. So far we do not have a unified 
framework which would allow for a systematic comparison on the different scales currently tested by cosmology. 



GR 



MAG 



Branework 



Bimetric 



SG 



Newtonian 



Scalar 
Tensor 



f(R) 



NGT 



Brans-Dicht' 



Gravity theories 



Cosmological 
concordance model 



BBN 



Weak 



(Parametrized) 
Post-Newtonian 
cosmological formulation 



Alternative cosmological 
models 



Cosmological tests 



CMS 



LSS 



SNIa 



FRUb 



Strong 



Sensing J. tensing 



Micro 
tensing 



X-ray 
cluster 



Rotation 
curves , 



FIGURE 2. The great advantage of a systematic approximative framework is depicted here for a future parametrized post- 
Newtonian approximation in cosmology. In contrast to the current strategy of testing different gravitational theories by working out 
single cosmological tests one could quickly asses the viability of a given theory by comparing its predictions for some cosmological 
PPN parameters. As described in the text, we need a unified PM and PN inspired framework to cover all of the cosmological tests. 
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we would be glad to work out the optimal form of the 
field equations and equations of motion for numerical 
treatment. 

PROBLEMS & OUTLOOK 

Let us summarize that currently we do not have a quan- 
titative estimate for the contribution of post-Newtonian 
effects to some very crucial cosmological tests. An im- 
portant example being the formation of structures in the 
universe. Hence it is of utmost importance to study these 
effects numerically in a simulation in order to assess up 
to which level we can trust current simulations which 
only make use of the Newtonian limit. Up to 1PN order 
we derived all of the necessary equations and would be 
glad to assist other groups in the numerical implementa- 
tion of these equations. 

On the theoretical side one of the open challenges of a 
cosmological PN approximation is linked to the need for 
conditions which clearly specify the realm of applicabil- 
ity of such an approximation. That means that we have 
to find a way to control the errors which we inevitably 
introduce when we treat a relativistic gravity theory in 
a Procrustean way. Another open issue in the cosmolog- 
ical case is the possible emergence of divergencies and 
thereby breakdown of the approximation at higher post- 
Newtonian orders. Such a behavior is already well known 
in the case of the PN approximation for isolated gravitat- 
ing systems [28]. Although it is unlikely that we need 
the improved accuracy of higher PN orders in any cos- 
mological test in the near future, a thorough theoretical 
treatment should also address this fundamental issue. 

As depicted in figure 2, the cosmological model 
building process would also greatly benefit from a 
parametrized post-Newtonian (PPN) framework. Such a 
scheme would allow for a rapid comparison of different 
theories by comparing their predictions for a set of stan- 
dardized cosmological PPN parameters. To cover all cos- 
mological tests one should aim for a framework which 
unifies PM and PN techniques. 

In view of a possible future PPN scheme for cos- 
mology the question emerges whether General Relativ- 
ity or, more precise, a metric theory of gravitation, is the 
appropriate foundation for such a scheme. As becomes 
clear from figure 2, there is no lack of alternative theo- 
ries which have been used to explain current cosmolog- 
ical observations. Some of these theories deviate signif- 
icantly from GR. Therefore, our choice of the structure 
underlying a future PPN inspired cosmological frame- 
work should be guided by the need to incorporate as 
many richer theories as possible. 
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